Can MSSM with light sbottom and light gluino survive Z-peak constraints ? 
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In the framework of minimal supersymmetric model we ex- 
amine the Z-peak constraints on the scenario of one light sbot- 
tom (2 ~ 5.5 GeV) and light gluino (12 ~ 16 GeV), which 
has been successfully used to explain the excess of bottom 
quark production in hadron collision. Such a scenario is found 
to be severely constrained by LEP Z-peak observables, espe- 
cially by Rb, due to the large effect of gluino-sbottom loops. 
To account for the Rb data in this scenario, the other mass 
eigenstate of sbottom, i.e., the heavier one, must be lighter 
than 125 (195) GeV at 2cr (3<r) level, which should have been 
produced in association with the lighter one at LEP II and 
will probobaly be within the reach of Tevatron Run 2. 
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Introduction Although the standard model (SM) has 
been successful phenomenologically, it is generally be- 
lieved to be an effective theory valid at the electroweak 
scale and some new physics must exist beyond the SM. 
This belief was seemingly corroborated by some exper- 
iments, such as the recent measurement of muon g — 2 
|jj and the evidence of neutrino oscillations Q. Among 
various speculations of new physics theories, the minimal 
supersymmetric model (MSSM) is arguably a promising 
candidate and has been intensively studied in the past 
decades. 

The non-observation of any sparticles from direct ex- 
perimental searches suggested heavy masses for sparticle 
spectrum. However, there have been a lot of analysis |^] 
which argue that a very light sbottom and light gluino 
(with mass of a few GeV) may have escaped from the 
direct experimental searches. It is intriguing that a light 
sbottom may require a light gluino, as analyzed in the 
last reference in J|. A recent analysis showed that 
a light sbottom (&i) with mass comparable with bottom 
quark is still allowed by electroweak precision data if its 
coupling to Z boson is small enough. A study by Berger 
et al || found that the scenario of MSSM with one light 
sbottom (2 ~ 5.5 GeV) and a light gluino (12 ~ 16 
GeV) can successfully provide an explanation for the 
long-standing puzzle that the measured cross section of 
bottom quark production at hadron collider exceeds the 
QCD prediction by about a factor of 2 ||. They also 
argued that such a scenario is consistent with all experi- 
mental constraints on the masses and couplings of spar- 
ticles. 



We note that the previous examinations j|] on Z-peak 
constraints focus on the direct production of a light sbot- 
tom followed by its decay similar to the bottom quark. 
Then by fine-tuning the mixing of left- and right-handed 
sbottoms, the coupling of Z-boson to the lighter mass 
eigenstate of sbottom (bi) can be sufficiently small so 
as to avoid the Z-peak constraints. It is noticeable that 
when sbottom b\ and gluino are both light, as was used to 
explain the excess of bottom quark production in hadron 
collision (5) , gluino-sbottom loops may cause large effects 
in Zbb coupling []. Therefore, in such a scenario, it is im- 
portant to reexamine the loop contributions to Zbb cou- 
pling and further, the Z-peak constraints. This is the aim 
of this letter. Through explicit calculations, we do find 
that gluino-sbottom loops comprising of sbottoms and a 
light gluino cause large effects on Z-peak observables. To 
account for the Rb data, subtle cancellation between b\ 
loops and 62 loops is needed, which can be realized by 
requiring the mass splitting between two sbottoms not 
to be too large. Numerical results show that for b\ with 
mass of 2 ~ 5.5 GeV, 62 must be lighter than 125 GeV 
and 195 GeV at 2a and 3<r level, respectively. 

Calculations We start the calculations by writing 
down the sbottom mass-square matrix [0 
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mass terms for left-handed squark doublets Q and right- 
handed down squarks, respectively. is the coefficient 
of the trilinear term H\QD in soft- breaking terms and 
tan /? = V2/V1 is the ratio of the vacuum expectation 
values of the two Higgs doublets. By diagonalizing the 
sbottom mass-square matrix, one obtains the physical 
mass eigenstates &i 2 
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1 Previous calculations of SUSY loop effects on Zbb coupling 
focused on rather heavy squarks and gluino and thus obtained 
very small effects M . 
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where 9 is the mixing angle of sbottoms. In our following 
analyses we take the sbottom masses and the mixing an- 
gel as free parameters since they are independent of each 
other and determined by SUSY parameters M~ , M? 

J r b L L b RR 

and M? . 

OLR 

The coupling of Z-boson to sbottoms is given by 



V^(Zb i b* j )=ieO i j(pi+p 2 y 



(3) 



and bj, respec- 
= v b + at, cos 29, 
O21 = —at sin 29. 



where Pi 2 are the momentum of b 
tively. Oij are defined as On = 
O22 = Vb — o-b cos 26 and 0\ 2 
Here Vb — 1/(4 sin 0vp cos 9\y)(1 — ^sin 2 9w) and a b = 
1/(4 sin 6w cos 6w) are the vector and axial vector cou- 
plings of Zbb, respectively. 

Apparently, a light sbottom b\ (a few GeV) can af- 
fect Z-peak observables in two ways: (1) the direct 
pair production of b\ through Zb\b\ coupling, as dis- 
cussed in 0]; (2) the loop effects of b\. If gluino is 
also light (12 ~ 16 GeV), then the loop effects are 
mainly from gluino-sbottom loops in Zbb vertex, which 
comprise a light gluino g and sbottoms, as shown in 
Fig. 1. It should be noted that even if the direct pair 
production of b\ is avoided by tuning the mixing angle 
I cos 6*| ~ ^/2/3sin Ow — 0.38 to set Zb\b\ coupling to be 
zero (On ~ 0), Zb\b* 2 and Zb 2 b\ couplings still exist and 
the irreducible loops shown in Fig. 1(b) make contribu- 
tions. It should also be noted that the self-energy loops 
in Fig. 1(a) involve only SUSY QCD interactions, i.e., 
gluino-sbottom-bottom couplings, which are not affected 
by the zero Zb±b* coupling. 
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FIG. 1. The gluino one- loop diagrams for Zbb. 

Using dimensional regulation and adopting the on-shcll 
renormalization scheme for the calculation of Fig. 1 , we 
obtain the effective Zbb vertex 
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lff {Zbb) = ie \-f^{v b - a b -f 5 ) + 7^ [Flip + ^27^75 

+i F 3 a^k" + 1 F 4 a^k»j 5 }} . (4) 

Here Fi are form factors originated from loop corrections, 
given by 

2 

Pi - 2 J2Oij{-Ar.m b rn- g (C (i,j) + 0^(1,3)) 
+A± [ml (C n (i,i) + C 21 (i,j)) + C 24 (i,j)] } 



F 2 = 
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+Vb5Z v + a b SZA, (5) 
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2 Y, Oii5±C 24 (i,i) - v b 6Z A - a b SZ v , (6) 



-J2 O t3 {A+m b [Cn(i,j)+C 2 i(i,j)] 

-m § Ar.[Co(i,j) + C u (i,j)]}, (7) 
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^4 = $30«{B+m 6 [2Cia(i,i)-Cii(i,j)-au(i,i) 
+2C 23 (i,j)]+m s B7-. [C Q {i,j) + C 11 {i,j)\] , (8) 
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Here B 0A (j) = Bo t i(-p b ,m g ,m b .) and C ,nm(hj) = 
Co tnm (—pb,k,mg,m b ,rni.), with p b and k denoting the 
four-momentum of b quark and Z boson respectively, are 
the Feynman loop integral functions and their expres- 
sions can be found in |9). Other constants appearing 
above are defined by Af- — diCij ± bibj, B^ = aibj ± ajbi, 
ai.2 = (sin 9 =F cos 6*) j\[2 and b\^ 2 = (cos 9 ± sin 0)/v / 2. 

Numerical results Let's now evaluate the effects of 
the above corrections to Z peak observables. We start 
with R b = T{Z -> bb)/T(Z -> hadrons). From Eq. (|) 
we obtain the contribution to R b 



SR b = R^ M (l~RDA SVSY , 
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with j3 = y/1 - Am b /m%. 

To obtain numerical results, we set input parameters as 
PI Rl xp = 0.21642 ±0.00065, R$ M = 0.21573 ± 0.0002, 



sirf 6 W = 0.2312, a s (m z ) = 0.1192, m z = 91.188 GeV 
and m b = 4.75 GeV. We will vary m g in the range 12 ~ 
16 GeV and in the range 2 ~ 5.5 GeV as was used 
in |. 

For my = 3.5 GeV and m g = 14 GeV, we present 5R b 
versus mg in Fig. ||. In addition to cos 9 = ±0.38 which 

leads to zero Zb\b\ coupling and hence avoids large rate 
of direct pair production of b\ [|J, we also plotted the 
curves for cos6> = ±0.30 and ±0.45. From the figure, one 
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sees that the contributions to Rb are negative in all the 
parameter space we have investigated. One can also see 
that the negative cos# gives larger contributions than 
positive one and as |cos#| increases, the contributions 
become more sizable. 
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FIG. 2. 5Rb as a function of m^ 2 for mi = 3.5 GeV and 
rrig = 14 GeV. The corresponding region above each hori- 
zontal line is allowed by LEP Rb data at 2a and 3<r level, 
respectively. 



Comparing with the experimental bounds shown in 
Fig. one learns that even in the favorable case of posi- 
tive cos 9, the contribution to Rb is too large to be allowed 
at 3 a level if mt > 200 GeV. Since the heavier sbottom 
has not been observed at LEPII, and it can in princi- 
ple be produced in association with the lighter one, its 
mass should probably be larger than about 200 GeV []. 
So, we conclude that the scenario of one light sbottom 
and light gluino faces severe challenge. As to the large- 
ness of the gluino-sbottom loop contributions, two main 
reasons may account for it. One is the large splitting be- 
tween mg i and , which leads to a weak cancellation 

between b\ and &2 contributions; the other is the lightness 
of sbottom b\ and gluino, which induces large self-energy 
contributions. To check our understanding, we fix 
and rrig but let approaches to mg 2 . Then we do find 
large cancellation occurs between different diagrams. 



We notice from Fig. || the intriguing feature that as 
increases, the effects get more sizable. This can be 

understood as the weaker cancellation between b\ and 62 

contributions when mz increases. To further understand 

2 

this behavior, we used approximate forms of B and C 
functions and found that in the limit to? 3> m? z > 

to? SR}, is roughly linear dependent on ln(mg2 /trip) 
and thus increases as mp / myi gets larger. Of course, this 
feature does not mean that SUSY QCD is non-decoupling 
from the SM. To check the decoupling property of SUSY 
QCD, we let all relevant sparticles (61, 62, g) become 
heavy and found the contributions drop quickly to zero. 
Actually, even for a light 61, SR b drops monotonously to 
zero when rrig get large, as shown in Fig. |^. 



CO 



2a bound , . . ^x****' 


3a bound . - ' 


: ,<^/ 




■'tM»' A/ 

: y 






cose= -0.3S 

cose= 0.38 


7 





FIG. 3. 



Mass of gluino (GeV) 

Same as Fig. ^[ but versus gluino mass for 



2 A detailed study may be needed to make this bound 
quantitative. 



m? = 3.5 GeV. 



Since in such a scenario with a light b\ of a few GeV, 
the &2 lighter than 200 GeV might be disfavored by LEP 
II experiment we fix = 200 GeV and cos 6 = 0.3 
and plot SRb versus in Fig. ^, where varies in 
the range 2 ~ 5.5 GeV and m g in 12 ~ 16 GeV, as used 
in Q to explain the excess of bottom quark production 
in hadron collision. We see that such a scenario is totally 
excluded by the LEP Rb data at 2cr level, while at 3tr 
level only a tiny corner with m g close to 16 GeV and 
close to 5 GeV is allowed. 

Let's next consider the effects on other Z-peak ob- 
servables: R Cl Re, Ab and A b FB . In our calculation of 
these observables, we neglect SUSY QCD correction to 
T(Z —>■ qq) (q ^ b) since the corresponding loops involve 
squarks q (q 7^ b) which arc assumed to be heavy Then 
the effects on all these observables stem only from the 
corrections to Zbb vertex in Eq.(|j). Since F\p are found 
to be much larger than F 3j 4, we neglect F 3i 4 in the cal- 
culation of Ab and A b FB . In Table 1, we show the effects 
on these observables including Rb- We see that gluino- 
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sbottom loop effects significantly enlarge the deviations 
of the predictions from the experimental values. 
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FIG. 4. Same as Fig. g, but versus mz for mg 2 = 200 
GeV. 



Table 1: Deviation of some Z-peak observables from exper- 
imental values. The MSSM predictions are obtained by in- 
cluding SUSY QCD contributions with = 3.5 GeV and 
m- g = 14 GeV. The SM predictions are taken from H. The 
values of mr are in units of GeV. 
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MSSM 


SM 


cos 6 = 0.30 


cos 9 = 0.45 


150 


200 


250 


150 


200 


250 


R b 


2.66(7 


3.22cr 


3.59cr 


4.49cr 


5.47cr 


6.16(7 


1.12(7 


R c 


-0.19cr 


-0.22cr 


-0.24cr 


-0.28cr 


-0.33(7 


-0.36cr 


-0.12cr 


Ri 


2.26ct 


2.66cr 


2.93cr 


3.65ct 


4.32cr 


4.83cr 


1.11(7 


A b 


-0.90cr 


-0.93cr 


-0.94cr 


-0.76(7 


-0.80(7 


-0.84a 


-0.64(7 


A b 

FB 


-3.25cr 


-3.28cr 


-3.30(7 


-3.06(7 


-3.12(7 


-3.16cr 


-2.90(7 



We should remind that in the calculation we only con- 
sidered the SUSY QCD loops, i.e., gluino-sbottom loops. 
Since we focused on a special scenario of the MSSM, in 
which there exist a very light sbottom (2 ~ 5.5 GeV) and 
a very light gluino (12 ~ 16 GeV), such gluino-sbottom 
loop effects are much larger than SUSY electroweak cor- 
rections Q. In fact, we re-calculated SUSY EW correc- 
tions to Rf, and found they are indeed small under the 
current experimental limits on the masses of charginos 
and stops. The dominant contributions from chargino 
loops are found to be positive (opposite to SUSY QCD 
corrections) and at the level of 10~ 4 , which are about one 
order smaller than our present SUSY QCD corrections. 

Conclusions From the above analyses we conclude 
that the scenario of the MSSM with one light sbottom 
(2 - 5.5 GeV) and light _gluino (12 - 16 GeV) can give 
rise to large effects on Zbb vertex through gluino-sbottom 



loops. Such effects significantly enlarge the deviations 
of some Z-peak observables, especially Rb, from their 
experimental data. To account for the Rb data in this 
scenario, the other mass eigenstate of sbottom, i.e., the 
heavier one, must be lighter than 125 (195) GeV at 2a 
(3ct) level, which should have been produced in associa- 
tion with the lighter one at LEP II and will probobaly 
be within the reach of Tevatron Run 2. 
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